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We study the effect of an external driving force on a simple stochastic reaction- 
diffusion system in one dimension. In our model each lattice site may be occupied 
by at most one particle. These particles hop with rates (1 ± f?)/2 to the right and 
left nearest neighbouring site resp. if this site is vacant and annihilate with rate 1 if 
it is occupied. We show that density fluctuations (i.e. the m th moments (N m ) of the 
density distribution at time t) do not depend on the spatial anisotropy ij induced by 
the driving field, irrespective of the initial condition. Furthermore we show that if 
one takes certain translationally invariant averages over initial states (e.g. random 
initial conditions) even local fluctuations do not depend on rj. In the scaling regime 
t ~ L 2 the effect of the driving can be completely absorbed in a Galilei transfor- 
mation. We compute the probability of finding a system of L sites in its stationary 
state at time t if it was fully occupied at time t$ = 0. 



PACS numbers: 05.40.+j, 05.50.+q, 02.50.Ga, 82.20.Mj 



1 Introduction 



Stochastic reaction-diffusion processes in one dimension have received a considerable 
amount of attention (for a vrief review see e.g. Despite their simplicity they 

show a very rich behaviour and some of the results obtained so far are of experimental 
relevance || . Another reason for their popularity is the wide range of applicability: 
Systems of this kind map to interface dynamics [[|, to polymers in random media 
@] or, quite close to every day life, to traffic problems ||. This list is by no means 
exhaustive, but provides already as it is ample motivation for the study of such 
models. Last, but not least, they map to well-known problems in many-body physics, 
particularly to integrable vertex models, and by now a considerable amount of exact 
and rigorous results have been obtained using this mapping 0. 

Here we consider a model defined on a ring of L sites with periodic boundary 
conditions where each lattice site may be occupied by at most one particle. These 
particles (denoted A) hop with rates (l±//)/2 to the right or left nearest neighbouring 
site resp. if this site is vacant (denoted 0) and annihilate with rate A if it is occupied: 

A0 -> U (1+W2 (LI) 

(DA -> A0 (1 - 77)/2 (1.2) 

AA — > 00 A (1.3) 

(In this paper we study only A = 1.) This model, which is closely related to zero 
temperaure Glauber dynamics [171 |TD[, has been studied by a number of authors over 
the past few years |10j - [|l4j . In the limit A = the model reduces to the well-known 



asymmetric exclusion process |To| whereas for A — ► oo, equivalent to the absence 
of any diffusion, the model becomes equivalent to random sequential adsorption of 
dimers [16|]. A discrete time version of the model was studied in ||17|| . Physically, the 



exclusion principle corresponds to a hard-core on-site repulsion and the asymmetry 
in the left and right hopping rates may be thought of as the result of a field driving 
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the particles in one preferred direction. The pair annihilation finally takes into a 
acount the possibility of an effective short-range attraction leading to unbreakable, 
inert pairs of particles. 

For A = the effect of the driving field has been well studied and it turns out to 
be rather drastic. As opposed to the driven model without exclusion where the drift 



may be absorbed in a simple lattice Galilei transformation [LJ| , here the asymmetry 
leads to the formation of shocks |19|]. Not many exact results are known for the 
time evolution of the system, but the mapping to the six-vertex model has been 



shown to be useful by using Bethe ansatz [T20], |2T[ and related methods |22| . Among 
other things the dynamical exponent appearing in the dynamical structure function 
of the system was found to be z = 3/2 |2(J rather than z = 2 in the undriven 
exclusion process or the non-interacting driven system without exclusion. For non- 
zero, but small annihilation rates A a recent study based on scaling arguments, 
mean-field approaches, random walk considerations and numerical results has given 
a very different scenario: in the presence of pair annihilation the effect of the driving 



appears to be very small In addition to that, it is known that for A = 1 and 

an initially full lattice neither the time-dependent density nor the two-point density 
correlation function depend on the driving HlOfl . These observations make a more 
detailed study of the system desirable. 

It is the aim of this paper to provide some exact results concerning the driving 
in the presence of the reaction. In Sec. 2 we define the model in terms of a master 
equation written in a quantum Hamiltonian formalism. We shall focus on the case 
A = 1 because for that particular choice the model may be described in terms of free 
fermions and treated rigourosly. Exact expressions for various correlation functions 
(in principle, for all correlation functions) become readily available. It turns out 
(Sec. 3) that fluctuations in the total density of particles do not depend on the 
driving field at all. The same is true for some other average values, to be specified 
below. Furthermore, if one takes certain translationally invariant averages over 
initial states, e.g. random initial conditions, arbitrary local density correlations show 
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no dependance on the asymmetry rj. These results generalize earlier findings p3| . 
The general solution of the master equation and some further results concerning the 
scaling regime are given in Sec. 4. In Sec. 5 we summarize and discuss our results. 



2 The master equation in quantum Hamiltonian 
formulation 

2.1 Definitions 

We define the process in terms of a master equation for the probability f(n;t) of 
finding, at time t, any configuration n of particles in the system of L sites with 
periodic boundary conditions. Here n = {n.1,712, . . . where rii — 0,1 and 1 < 
i < L labels the sites of the lattice. An alternative possibility is to give the set 
{x\, X2, ■ ■ ■ ,xn} of occupied lattice sites. In this notation, the empty set represents 
the empty lattice and 1 < iV < L is the total number of particles in the configuration. 
We shall express the time evolution given by the master equation in terms of a 



quantum Hamiltonian H |24| . This is discussed in detail in a number of publications 
(for consistency of notation see e.g. ||) and we shall repeat only the essential 
elements of the mapping. The advantage of this approach is that there are standard 
methods of dealing with the resulting time evolution operator H. The applicability 
of these techniques, in the case at hand essentially a Jordan Wigner transformation 
and the representation of states in terms of a fermionic Fock space, does not arise 
naturally and obviously if the master equation is written down in standard form. 

The idea is to represent each of the 2 L possible configurations in X = {0, 1} L 
by a vector \n) (or \xi, . . . ,x^), with |0) = | ) being the empty state). The 
probability distribution is then mapped to a state vector 

!/(*)>= £/(ffi*)la> • (2.1) 

The vectors | n ) together with the transposed vectors ( n \ form an orthonormal basis 
of (C 2 )® L and the time evolution is defined in terms of a linear 'Hamilton' operator 
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H acting on this space of dimension 2 L 

^\f(t)) = -H\f(t)) . (2.2) 

A state at time t = to + T is therefore given in terms of an initial state at time to by 

|/(to + r> = e- HT |/(t > • (2.3) 

From ( |2.1| ) and ( |2.2|) and using /(n;t) = (n\f(t)) the master equation takes the 
form 

9 f(n;t) = -(n\H\f(t)) . (2.4) 



Of 

Note that 



where 



(s\f(t))=J2f(n;t) = l (2.5) 

nex 



£<n| (2-6) 



which expresses conservation of probability and which implies ( s \H — 0. 

Expectation values ( Q ) are calculated as matrix elements of suitably chosen 
operators Q. A complete set of observables are the occupation numbers = 0,1. 
Defining projection operators on states with a particle on site k of the chain as 

1 ( \ 

n» = 5 (l-^)=^ (2-7) 

one finds that the average density of particles at site k is given by ( rik ) = 
( s \nk\ f(t) ). Correlation functions ( • • • ) , i.e., the probabilities of finding 
particles on the set of sites {ki, . . . , kj}, are computed analogously. 

For later convenience we also introduce the operators = (af ± ia y k )/2. In our 
convention 

/ o\ 

{21 





1 



k 

creates a particle at site k when acting to the right, while 



4 = | H I (2.9) 



annihilates a particle at site k. Note that 

(s\s% = (s\n k and ( s \s% = ( s |(1 — n k ) . ( 2 -10) 

Introducing the ladder operators = J3fe =1 s k one may write 

(s| = (0|e 5+ . (2.11) 

Using the commutation relations for the Pauli matrices then yields (|2.10| ). 
Now we can define the process in terms of the quantum Hamiltonian 

L 

(2.12) 

k=l 

with the nearest neighbour reaction matrices 

1 + Ti t , _ \ 1 — 7] / , \ 

= [M 1 - n k+i) - sts k+1 J + — ^— (Jl - n k )n k+1 - s k sl +1 J 

+A (n k n k+1 - (2.13) 
= — ^— K - sit — ^— K+i - s fc s£ +1 ) 



Wfe 



This, together with ( |2.4j ) defines the process. The Hamiltonian may be written 
H = H s +r]H d where the driving part H d is given by H d = l/2J2 k =i(s k sl +1 -s k s k+1 ) 
and H s is the Hamiltonian for the system without driving. 

2.2 Fermion representation of H 

For A = 1 the Hamiltonian becomes bilinear in the creation and annihilation oper- 
ators sf. This suggests rewriting H by introducing fermionic operators through a 



Jordan- Wigner transformation [25]. We define 



Qk 


k 

= IK 


(2.14) 




i=l 




4 


= s k Q k -i 


(2.15) 


«k 


= Q fc _ a s fc 


(2.16) 
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satisfying the anticommutation relations {a k ,ai} = {4,4} = and {a\,a{\ = 8 k ,i- 
Note that because of the periodic boundary conditions for the Pauli matrices one has 
4+i = a \Qh and cll+i = QlO-i- Ql may be written Ql = {— 1) N where N = n k is 
the number operator. Since by the action of H the particle number changes only in 
units of two, Ql commutes with H and splits it into a sector with an even number of 
particles (Ql = +1) and into a sector with an odd number of particles (Ql = — 1). 
In terms of the fermionic operators one finds 

SfcSfc+i = a\ +l a k Sfc4+i = a\a k+1 (2.17) 
4 4+i = Ofc+iOfc SfeSfc+i = 44+i (2.18) 
n k = a\a k (2.19) 

and we arrive at the following expressions for H s and H d : 

1 L 

Hs = ~o {4+i°fc + alofc+i + 2a k+1 a k - 2n fe | 
1 k=i 

L 

-(A - 1) {afc+iafc - n k n k+1 } (2.20) 
fc=i 



Ha 



1 

k=l 

We also note that H s and i/^ commute if A = 1. Independently of A and rj the 
non-degenerate ground state of H, corresponding to the steady state of the system, 
is the totally empty state |0) in the sector with an even number of particles and the 
state \0) odd = 1/LEti \ k ) where 

one particle may be found with equal probability 
1/L anywhere in the lattice in the sector with an odd number of particles. 

Since expectation values (n k (t)) of the stochastic variables n k = 0,1 are 
given by the matrix elements (sln^exp (—Ht)\f) of the operator n k and since also 
(s\ exp (Hi) = (s\, one may introduce time-dependent operators 



0(t) = exp (Ht)Oexp -(Ht) 
6 



(2.22) 



and write (0(t)) = (s\0(t)\f) for an arbitrary operator O. According to the defini- 
tion ( p.22|) one has 

jO = 6= [H, O) (2.23) 

Among the quantities of interest are the density of particles in the system and 
fluctuations in this quantity. The m th moment of the particle number distribution 
is the expectation value (N m (t)) where N is the number operator N = Y^k=i a \ a k- 
Any O may be written as a product of the fermionic annihilation and and creation 
operators and it is therefore sufficient to study the time evolution of these operators. 
It is useful to introduce the Fourier transforms 

~i- L 

b p = "7fE e ~^ ( 2 - 24 ) 
v L fc=l 

i— L 

b l = -7fT, e ~~^4 (2-25) 

V L fc=1 

satisfying {b p ,b q } = {6j, 6j} = and {6j, 6J = V Inverting (p4|), (pD yields 

= ^rE e ^ & p (2-26) 
V L p 

„— in/4 „ ., 

4 = -T^E^J (2 . 27) 
V i/ p 

Thus the representation of the number operator in Fourier space is 

N = Y.b% (2.28) 

p 

Here the sum runs over all integers p — 0, . . . , L— 1 in the sector with an odd number 
particles and over the half odd integers p = 1/2, 3/2. . .L — 1/2 in the even sector. 
The Hamiltonian for A = 1 reads 

H s = ^|(l-cos^)fetfo p + sin^6_,6 p } (2.29) 

# d = -z]Tsin(^)& P & p (2.30) 
p L 

H s was already obtained in E3 



2.3 Initial states in the fermion representation 

Having discussed the representation of operators in terms of the fermionic operators 
we turn to the representation of initial configurations or states. In what follows we 
shall use the term '(initial) configuration' for a vector \ki, . . . , k^) or |n), i.e., for a 
simultaneous eigenstate of all n^. (Initial) states |/) = J2 n f(l±)\ll) are vectors which 
may be superpositions of such configurations, and normalized such that (s\f) = 1. 
The vector |0) is the vacuum state with respect to the annihilation operators, afc|0) = 
0. In spin language this is the ferromagnetic state with all spins up. Acting with 
creation operators yields 

4r--4j°) = \h,---,k N ) (h < k 2 < . . . < k N ) (2.31) 

A general translationally invariant iV-particle state is obtained by acting with 
products tf '"bp N on |0) where J2iPi = 0. A special class are those build by 
polynomials in the bilinear expressions B\ = &t p &t where p= 1/2, 3/2,..., (L-l)/2. 
Among these particular translationally invariant states are uncorrelated random 
initial conditions with an even number of particles. The uncorrelated random initial 
state with density p is the product measure 

\P) = ( 1 T P = E - P) L - N \N) (2.32) 

\ r J N=0 

where \N) is the iV-particle state where each configuration appears with equal 
weight. By projection on the sectors with even and odd particle numbers we obtain 

lpr en(o dd ) = i ^ QL 2p)L \p) (2.33) 

where we have used Ql defined in ( p.!4| ). 

From the representation (|2.29|) , (|2.30|) of H and using (B p ) 2 = (B^) 2 = one 



finds for the left zero energy eigenvector (s\ of H 

(s\ = (O|n(l + cot(^)6 p 6_ p ) +(0|6oll( 1 + cot (x) fop ' 6 - p ') 



S 



( 0| e ^( 1+cot (^KM + (0|6 e Ep '( 1+COt ( ¥ ) V6 - p ') 

^jeuen _|_ ^|odd 



[2M) 



E-t(f)4 = E (2.35) 



The product and sum resp. over p,p' run over p = 1/2,3/2, . . . , (L — l)/2 (even 
sector), and over p' — 1,2, . . . ,L/2 — 1 (odd sector). This in turn implies 

v 

iV! 
Hence 

|2JV> = Mfe cot (f) B ') |0> (2 ' 36) 

W "" = l + ( l-^ n((l-rt' + ^cot(g)^)|0> (2.37) 
The completely full lattice is simply given by 

|L) = n^|0) • (2.38) 

p 

In what follows we shall focus on the even sector. For a study of correlation 
functions (Sec. 4) it is useful to note 



even 



&t+cat(^V p ) =0 (2.39) 



which may be verified using the momentum space representation ( |2.34j ) of (s| 



3 Dynamics of the system on a translationally in- 
variant subspace 

In the last section we have introduced the subspace V generated by the operators 
BJj. Together with B p = b p b^ p they satisfy the algebra of the Pauli matrices, with 
[Bp, BU = b-pb^p — bpbp = 2C P playing the role of the a p matrix. I p = b p b p — b^ p $_ p 
commutes with all B^, B p , C p and acts as unit operator on this subspace. It satisfies 
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IpXp = Xplp = Xp for X p = Bp, Bjj, C p . These relation are easy to verify by using 
the anticommutation relations for bp and b p . 

This subspace is of interest for three reasons. Firstly, H can be written in terms 
of these operators: 

H = E ( x - cos ( ) Np ~ 2 sin (t) Bp ~ 171 sin (t) (/p " 1} (3 ' 40) 

Where we have introduced the number operator N p — 1 — 2C P . Therefore V is an 
invariant subspace of H. Secondly, a physically important class of initial conditions, 
namely random initial conditions, including the steady state and the fully occupied 
lattice, are in this subspace. Finally, some physically important expectation values 
are given by operators constructed from and B p . In particular, using ( |2.2£| ) one 
gets 

iV = £iV p (3.41) 

p 

Without further calculation we can now state the following results: 

1 ) The subspace V of dimension 2 L I 2 generated by B^ acting on the vacuum state 
|0) is an invariant subspace of H . On this subspace = 0, i.e., the driving has no 
effect on any correlation function if the system is at time t = in an initial state 
which is contained in V (e.g. random initial conditions). 

That this correct can be seen by observing that I p is the unit operator on this 
subspace which gives H d = 0. As a result, the state at time t does not depend on r\ 
which in turns implies that no correlation in that state can depend on the driving. 

2) The time evolution of operators O build from operators B^ and B p (e.g. the den- 
sity operator) does not depend on the driving, irrespective of the initial condition. 

This is again obvious since commutes with any such operator. Applying this 
result to powers N m of iV we find 
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3) The probability P(N;t) of finding precisely N particles at time t in the system 
does not, for any initial condition, depend on the driving. 



The same applies then of course also for the moments (N m (t)) of this distribution. 



For m = 1 this was shown in [10]. An interesting special case is Pf(Q;t) = (0\f(t)) 
giving the probability that the system has reached the steady state at time t from 
an initial state |/). To obtain this quantity we note that H ( |3.40|) restricted to the 



subspace V is a sum of 2 x 2 matrices h 



v 



( -sin27r»/L \ 

H = 2 ^[ n 1 9 ,t) (3 ' 42) 
P \0 1 - cos 2vrp/L J p 

Computing exp —Ht = Ilp ex P —h p t and taking the product of the matrix elements 
(0\B^(t)\0) then gives for random initial conditions (|2.37|) with an even number of 
particles the following exact expression: 



For an initially full lattice this simplifies to 



e -2t(l-cos (2tt P /L))\ 



(3.43) 



P^Q-t) = J] (l - e -M0.-«*(2*P/L))j I^po e -L/V^t (3 44) 

V 

4 Full solution of the master equation 

The master equation is solved if one has found explicit expressions for a\{t) and 
afc(t) or their Fourier transforms. The knowledge of these quantities allows then the 
explicit calculation of any correlation function for any initial condition, including 
multi-time correlators giving arbitrary conditional probabilities. We will compute 
a\(t) and a^it) and study the effect of the driving on these quantities for late times. 
Applications to specific correlation functions of interest will be given elsewhere. For 
definiteness we study only the sector with an even number of particles and we also 
assume, as in the preceding section, L even. 
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It is interesting to study first the differential equation satisfied by the local 
density (n k (t)). Differentiating with respect to time one finds 

J t (Mt)) = ~ (K+iW) + (n k -i(t)) - 2(n k (t))) - | «n fc+ i(t)> - (n k ^(t))) 

-(A + v)(nk-i(t)n k (t)) - (A - r])(n k (t)n k+1 (t)} (4.45) 

In the linear terms one recognizes a lattice Laplacian and lattice derivative respec- 
tively. For A = the nonlinear term is also a lattice derivative and the equation is a 



discrete form of Burgers equation |26] describing the evolution of shocks for rj 7^ 0. 
Any A > will result in a strong dampening of the amplitude and the question 
arises to which extent the nonlinear effects associated with the driving continue to 
play a role. For A = 1 this question was partially answered in the last section, the 
result being, somewhat surprisingly, that, for the initial conditions considered, there 
are now such effects at all. Here we study the form of local correlations for arbitrary 
initial conditions. 

The simple form of H in Fourier space suggests studying &I(t) and b p (t) rather 
than a k (t) and a k (t). Applying (|2.23|) gives a set of two coupled ordinary differential 
equations^ 



ffiit) = e p &t(t)+2sin(^)Mt) (4-46) 
d b p (t) = -e p b p (t) (4.47) 



dt 



solved by 



bl{t) = (bl + cot ^ (l - e -^ +e -^) b^ (4.48) 
bJt) = e- e *X (4.49) 



with 6j(0) = 6j, 6p(0) = b p and 



2vrp . . 2np 
e p = 1 — cos — irj sin — - — (4.50) 



1 Because of the boundary conditions these equations describe the time evolution of the creation 
and annihilation operators only when applied to products with an even number of operators. This 
is not a restriction as all expectation values (n kl ■ ■ ■ n^ N ) are of this form. 
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From ( [4.48|) and ( |4.49| ) one obtains 



n k (t) = jT,e 2mk{p - p ' )/L bt(t)bA^ (4.51) 
p,p' 

and therefore an explicit expression of any correlation function at time t in terms of 



correlators in the initial state. In particular, one may use ( |2.39|) to obtain 



(Mt)} = t E e^-^/^^+v)* cot ^(V&p) (4-52) 

P,P> 

This solves the initial value problem for the average density. 

Returning to arbitrary time-dependent correlators we note that Q4.48|) - ( f4.50|) 



demonstrate the impact of the driving on the system in the scaling regime. For 
times t ^> L 2 all correlations decay exponentially with correlation time r = 2L 2 /tt 2 
resulting from slowest mode p = 1/2. For times in the scaling regime t ~ L 2 and L 
large, one may approximate e p by 

2ttz 2tt 2 „ , . 

€p ~ ~^~zr p + jtp ( 4 - 53 ) 

In other words, the effect of the driving may be absorbed in a Galilei transformation 
r i r i + where = k{/L are the scaled space coordinates appearing in the 
correlation function. Apart from that there are no other effects, unlike in the absence 
of annihilation where the driving induces the evolution of shocks. For arbitrary 
translationally invariant initial conditions the dependence of correlation functions 
on 7] vanishes completely in the scaling limit. 

5 Conclusions 

We have studied the effect of driving in a simple reaction-diffusion system where 
exclusion particles hop with rates (1 ± ^)/2 to the right and left respectively if the 
nearest neighbour sites are empty and which are annihilated in pairs with rate 1 if 
the nearest neighbour site is occupied. We obtained the following results: 
(1) Certain translationally invariant time-dependent expectation values including 
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the m th moments (N m ) of the particle number distribution do not depend on the 
driving parameter rj, regardless of the initial condition. 

(2) Arbitrary time-dependent expectation values of occupation numbers (n^ ■ ■ ■ Uk N ) 
do not depend on rj if at time t = one takes certain translationally invariant aver- 
ages over initial states, e.g. uncorrelated random initial conditions. 

(3) In the scaling regime t w L 2 the effect of the driving can be completely absorbed 
in a Galilei transformation. 

(4) The probability that the system has reached its steady state (the empty lattice) 
at time t from an uncorrelated random initial condition with density p and an even 
number of particles is given by ( |3.43|) . 

These results have been obtained for a annihilation rate A = 1. It would be 
very interesting to study the system for other values of A in order to understand 
the transition to the limiting cases A = and A = oo. For A = (the asymmetric 
exclusion process) the introduction of driving has a very strong effect, it causes the 
evolution of shocks from local inhomogeneities and it is not clear to which extent 
these effects survive in the presence of annihilation. Our calculation shows that 
for A = 1 they are completely absent in the scaling regime. It is interesting to 
note that it is not the exclusion principle as such which is responsible for the non- 
linear behaviour of the asymmetric exclusion process, but the strength of the pair 
interaction between neighbouring particles. In stochastic language this is the rate of 
change of a pair of neighbouring particles A compared to the diffusion rates (l±r/)/2. 
A satisfactory understanding of this interplay remains an open problem. Intuitively 
one would expect the nonlinear behaviour to vanish for late times for any A / 
as the system will then be almost empty and therefore effectively non-interacting. 
This is supported by various arguments put forward in Ref. |nj and by the known 
dynamical exponent z = 2 of the asymmetric exclusion process in the regime of 
low (infinitesimal) densities |L2] (as opposed to z = 3/2 for finite densities). But 
a qualification of the expression 'late times', and, more importantly, what happens 
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before such late times remains open to debate. 

Another interesting open problem is the approach to local equilibrium, i.e. to a 
local extended region of empty sites. The first question to be asked in this context is 
the probability of reaching a state with, say, M empty neighbouring sites. The next 
question then is how this state further evolves in time. It is of course not stationary 
as particles are injected and absorbed by diffusion at the boundaries of this region, 
but with time-dependent (vanishing) rates. For this problem the driving will make 
a difference as the injection of particles at the left boundary of this region will be 
stronger than the loss of particles (if particles move preferredly to the left). At the 
right boundary of this region the situation will be reversed, the loss will exceed the 
gain. The result of Sec. 4 show that the problem with a drift may be obtained from 
the symmetric problem by a Galilei transformation and therefore gives a partial 
answer to this question. Eq. ( [3.431 ) gives the probability of finding the system in 
global equilibrium which does not depend on the driving. 

6 Acknowledgments 

This work was supported by a grant of the European Community under the Human 
Capital and Mobility program. 

References 

[1] V. Privman, Dynamics of N onequilibrium Processes: Surface Adsorption, 
Reaction-Diffusion Kinetics, Ordering and Phase Separation preprint |cond-mat 
93T2079| 

[2] Mac Donald, Gibbs, Pipler, Biopolymers 6, 1 (1968); R. Kopelman, C.S. Li and 
Z.-Y. Shi, J. Luminescence 45, 40 (1990); R. Kroon, H. Fleurent and R. Sprik, 
Phys. Rev. E 47, 2462 (1993). 



15 



[3] J. Krug and H. Spohn, in: Solids far from Equilibrium, ed. C. Godreche, Cam- 
bridge, Cambridge University Press (1991) and references therein. 

[4] M. Kardar and Y. C. Zhang, Phys. Rev. Lett. 58, 2087 (1987). 

[5] K. Nagel and M. Schreckenberg, J. Phys. I, 2, 2221 (1992); A. Schadschneider 
and M. Schreckenberg, J. Phys. A 26, L679 (1993). 

[6] I. Peschel. V. Rittenberg and U. Schultze, Nucl. Phys. B 340, 655 (1994); G. 
Schiitz, J. Stat. Phys. (in press); and references therein. 

[7] R. J. Glauber, J. Math. Phys. 4, 294 (1963). 

[8] A. A. Lushnikov, Phys. Lett. A 120, 135 (1987). 

[9] J. L. Spouge, Phys. Rev. Lett. 60, 871 (1988). 

[10] M. D. Grynberg, T. J. Newman and R. B. Stinchcombe, Phys. Rev. E 50, 
957 (1994); M. D. Grynberg and R.B. Stinchcombe, Phys. Rev. Lett. 74, 1242 
(1995). 

[11] J. G. Amar and F. Family, Phys. Rev. A 41, 3258 (1990); F. Family and J. G. 
Amar, J. Stat. Phys. 65, 1235 (1991). 

[12] M. Henkel and G. M. Schiitz, Physica A 206, 187 (1994). 

[13] F. C. Alcaraz, M. Droz, M. Henkel and V. Rittenberg, Ann. Phys. (New York) 
230 . 250 (1994). 

[14] V. Privman, E. Burgos and M. D. Grynberg, Multiparticle Reactions with 
Spatial Anisotropy, |cond-mat 9410103] . 

[15] Zia and B. Schmittmann, Statistical mechanics of driven diffusive system, to 
appear in Phase Transitions and Critical Phenomena, eds. C Domb and J. 
Lebowitz (Academic, London). 



16 



[16] J. W. Evans, Rev. Mod. Phys. 65, 1281 (1993) 
[17] V. Privman, J. Stat. Phys. 72, 845 (1993). 

[18] M. Henkel and G. M. Schiitz, Int. J. Mod. Phys. B 8, 3487 (1994). 

[19] B. Derrida, S. A. Janowsky, J. L. Lebowitz and E. R. Speer, Europhys. Lett. 
22, 651 (1993). 

[20] L.-H. Gwa and H. Spohn, Phys. Rev. Lett., 68, 725 (1992), Phys. Rev. A, 46, 
844 (1992). 

[21] G. M. Schiitz, J. Stat. Phys. 71, 471 (1993). 

[22] G. M. Schiitz, Phys. Rev. E 47, 4265 (1993); S. Sandow and G. M. Schiitz, 
Europhys. Lett. 26, 7 (1994). 

[23] G. M. Schiitz, Oxford preprint. 

[24] L. P. Kadanoff and J. Swift, Phys. Rev 165, 310 (1968); M. Doi, J. Phys. A 9, 
1465, 1479 (1976); P. Grassberger and M. Scheunert, Fortschr. Phys. 28, 547 
(1980). 

[25] P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 

[26] J. M. Burgers, The non-linear diffusion equation, Riedel, Boston (1974). 



17 



